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Abstract 

The H atoms inside minihalos (i.e. halos with virial temperatures T v ; r < 10 4 K, in 
the mass range roughly from 1O 4 M0 to 10 8 M Q ) during the cosmic dark ages in 
a ACDM universe produce a redshifted background of collisionally-pumped 21-cm 
line radiation which can be seen in emission relative to the cosmic microwave back- 
ground (CMB). Previously, we used semi-analytical calculations of the 21-cm signal 
from individual halos of different mass and redshift and the evolving mass function 
of minihalos to predict the mean brightness temperature of this 21-cm background 
and its angular fluctuations. Here we use high-resolution cosmological N-body and 
hydrodynamic simulations of structure formation at high redshift (z > 8) to com- 
pute the mean brightness temperature of this background from both minihalos and 
the intergalactic medium (IGM) prior to the onset of Lya radiative pumping. We 
find that the 21-cm signal from gas in collapsed, virialized minihalos dominates over 
that from the diffuse shocked gas in the IGM. 
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1 Introduction 



Neutral hydrogen atoms in the early universe can be detected in absorption or 
emission against the cosmic microwave background (CMB) at redshifted radio 
wavelength 21 cm, depending upon whether their spin temperature Ts is less 
than or greater than that of the CMB, respectively. Minihalos which form 
during the " dark ages" have densit y and temperature high enough to appear 
m emission fjlliev et al.ll2"n02l . l2003h . The intergalactic medium (IGM), on the 



other hand, appears in either emission or absorption. New radio telescopes are 
being designed to detect this 21 cm signal. 



Iliev et al. ( 20021 ISFM hereafter) showed that atomic collisions inside mini- 



halos are sufficient to decouple the spin temperature of the neutral hydrogen 
from the CMB temperature. They predicted the mean and angular fluctuations 
of the corresponding 21cm signal by a semi- analytical calculation based upon 
integrating the individual minihalo contributions for different halo masses and 
redshifts, over the evolving statistical distribution of minihalo masses in the 
ACDM universe. The fluctuations, the y found, are substantial en ough to be de- 
tectable with future radio telescopes. Furlanetto fc Loebl ( 20041 ). on the other 



hand, have recently suggested that the gas in the diffuse IGM (prior to the 
onset of Lya radiative pumping) is also capable of producing the 21cm emis- 
sion signal and that the IGM contribution to the mean signal will dominate 
over that from gas inside minihalos. 

In order to quantify these effects, we have computed the 21 cm signal both from 
minihalos and the IGM at z > 8, using high-resolution cosmological N-body 
and hydrodynamic simulations of structure formation. We use a flat, ACDM 
cosmology with matter density parameter Q m = 0.27, cosmological constant 
= 0.73, baryon density f2& = 0.043, Hubble constant H = 70kms _1 Mpc _1 , 
ash- 1 = 0.9 and the Harrison-Zel'dovich primordial power spectrum. 



2 THE CALCULATION 



2. 1 Basics of 21 cm Radiation Background from the Dark Ages 



Foreground emission or absorption by neutral hydrogen atoms at redshift z is 
seen against the CMB at redshifted wavelength, 21(1+2;) cm. The spin temper- 
ature (Ts) of a hydrogen atom determines whether the signal is in emission or 
absorption. Emission occurs when Ts > Tcmb, while absorption occurs when 
Ts < Tcmb- Ts can deviate from Tcmb in various ways. A hydrogen atom 
can 1) absorb a 21 cm photon from CMB (CMB pumping), 2) collide with 
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another atom (collisional pumping) and 3) absorb a Lya photon to make a 
Lya transition, then decay to one of hyperfine 21 cm levels (Lya pumping). 

During the dark ages, when there is no Lya pumping, the mean 21 cm signal 
against the CMB will be zero at z > 150, then in absorption at 20 < z < 
150 mostly due to unperturbed IGM, and finally in emission at z < 20 due 
primarily to minihalos. We restrict ourselves to regions such that the Lya 
pumping is negligible even after sources turn on at z < 20. In other words, we 
only consider collisional pumping in minihalos and the IGM at 8 < z < 100. 



2.2 Semi-analytic calculations 



Here we briefly summarize the semi- analytical calculation of the 21 cm signal 
from minihalos by ISFM and from the unperturbed IGM. For minihalos, the 
differential brightness temperature, 

° T b = — — ( 1-e )' ^' 



is averaged over individual minihalos to give the mean differential antenna 
temperature 5T b (for detail, see ISFM). On the other hand, the unperturbed 
IGM of the universe has a kinetic temperature smaller than the CMB temper- 
ature at z < 100, resulting in an absorption signal until collisional pumping 
becomes negligible at z ~ 20 (Fig. [1]) . 



2.3 Numerical Simulations 



We have run a high resolution cosmological N-body and gasdynamic simula- 
tion to derive the effects of gravitational collapse and hydrodynamics on the 
predicted 21 cm signal from the high redshift universe. Our computational 
box has a comoving size of 0.7 Mpc, with 1024 3 cells and 512 3 dark matter 
particles, which is optimal for adequately resolving both the minihalos and 
small-scale struct ure-formation shocks. We have used the code described in 
Rvu et al. which uses the particle-mesh (PM) scheme for calculat- 



ing the gravity evolution and an Eulerian total variation diminishing (TVD) 
scheme for hydrodynamics. 

In order to calculate the minihalo and the IGM contributions to the total 
differential brightness temperature, one needs first to identify the halos in the 
simulation volume. We ident ified the halos using a friends-of-friends (FOF) 
algorithm ( Davis et al. 1985f ) with a linking length parameter of b = 0.25. 



The rest of the signal is considered to be from the IGM. 
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Fig. 1. Analytical prediction for the mean 21 cm differential brightness temperature 
due to collisionally-decoupled minihalos and an unperturbed IGM. Shown are the 
results based on the Press-Schechter (solid) and the Sheth-Tormen (dotted) mass 
functions for halos and the contribution from the IGM gas with cosmic mean density 
and temperature (dashed). In the bottom panel we show the minihalo collapsed 
fraction, again based on the Press-Schechter (solid) and the Sheth-Tormen (dotted) 
mass functions. 

In Figure 2, we show projected maps for three different redshifts. Each con- 
tribution from minihalos and the IGM is also plotted. It clearly shows that 
the signal is in absorption due to the IGM at z > 20, and after crossing the 
moment when a mixture of emission from minihalos and absorption from the 
IGM exists (z ~ 20), the signal is predominantly in emission both from mini- 
halos and the IGM at z < 20. In Figure 3, we show the total signal as well as 
minihalos and the IGM contributions. 



2.4 Semi- Analytical Calculation of the Halo Contribution 

Our numerical simulations have sufficiently high resolution to find all of the 
minihalos in the computational box, as well as the large-scale structure forma- 
tion shocks, but not to resolve the internal structure of the minihalos them- 
selves. As ISFM have shown, in order to obtain the correct 21-cm signal from 
minihalos one needs to do a full radiative transfer calculation through each 
individual minihalo density profile since, unlike the IGM gas, minihalos have 
a non-negligible optical depth at the 21-cm line. Hence, we can refine our es- 
timate of the minihalo contribution to the total 21-cm signal by combining 
our numerical halo catalogues with the semi-analytic calculation of individual 
minihalo contribution as found by ISFM. 

In their approach, the gas density of each minihalo is assumed to follow the 
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Fig. 2. Map of the 21 cm signal obtained from our high resolution simulation. Rows, 
top to bottom, show redshifts z=30, 20, and 10. Columns, left to right, represent 
contributions from minihalos, the IGM and the total signal, respectively. Note that 
the scale is linear for the upper two rows of images, but logarithmic for the bottom 
row. 



TIS profile, of Iliev fc Shapiro! ( 2001 ). the radiative transfer calculation is per- 
formed for different impact parameters, and then finally the face-averaged 
is calculated (see ISFM for details), according to 



jr. c(1 + z) 



Mn 



J Au ee 5T biVQ A-^dM. 



(2) 



The halo mass function, dn/dM, is provided here, not by the analytical PS 
mass function as in ISFM, however, but by the halo catologue we construct 
from the simulation. Each individual halo contribution, Au e s5Tb t VQ A, depends 
on its mass and redshift of formation (ISFM). Once we calculate Az/ eff <5T bj! , A, 
we then obtain the halo contribution using equation (2). 
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Fig. 3. 21 cm mean brightness temperature evolution. (a)(left) Evolution of the 
total 21-cm signal vs. redshift. All data points are directly calculated from our 
simulation box, with the assumption that optical depth is negligible throughout the 
box. (b) (right) ST^ vs. redshift below z = 20. Plotted are the contributions from 
minihalos (labeled Halos, circles), the IGM (triangles) and the total (squares). 
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Fig. 4. Semi-analytical minihalo signal vs. IGM signal. The 21 cm flux from each 
halo in the simulation is found by modeling the halos as described in § 2.4, to 
estimate the halo 21 cm signal more accurately. Same notation as in Figure 3. The 
21 cm minihalo emission increases compared to the raw minihalo signal in figure 3. 
The IGM signal remains the same. 



3 Minihalos vs. the IGM 



We find that the emission signal at z < 20 is dominated by minihalos. For 
8 < z < 14, about 70 — 80% of the total signal is contributed by minihalos. 
The minihalo contribution peaks at z ~ 16 — 18, because the IGM is in both 



6 



emission and absorption. At z ~ 20, the IGM is mostly in absorption, with an 
absolute amplitude which is comparable to the minihalo emission signal. Only 
when the collapsed fraction of baryons into minihalos is negligible (z > 20), 
does the IGM make the dominant contribution to the total signal, which is in 
absorption. Note that the total signal at z > 20 is slightly stronger than the 
analytical value for the unperturbed IGM, because clumping of the IGM makes 
the collisional coupling of the spin temperature to the kinetic temperature 
slightly stronger. 
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Fig. 5. Fractional contributions of minihalos and diffuse IGM gas to the total 21-cm 
background. The top panel shows the results obtained directly from simulations (CI: 
triangle, long-dashed; C2: square, short-dashed; C3: pentagon, dotted; C4: circle, 
solid). The bottom panel shows the results which were semi-analytically refined 
(§ 2.4; point- and line-types follow those of the top panel). 



4 Conclusions 

We have performed a cosmological N-body and hydrodynamic simulation of 
the evolution of dark matter and baryonic gas in the ACDM universe at high 
redshifts (8 < z < 100). With the assumption that radiative feedback ef- 
fects from the first light sources are negligible, we calculated the 21 cm mean 
differential brightness temperature signal. The mean global signal is in absorp- 
tion against the CMB above z ~ 20 and in overall emission below z ~ 18. At 
z > 20, the density fluctuations of the IGM gas are largely in the linear regime, 
and their absorption signal is well approximated by the one that results from 
assuming uniform gas at the mean adiabatically-cooled IGM temperature. At 
z < 20, nonlinear structures become common, both minihalos and clumpy, 
warm, mildly nonlinear IGM, resulting in an overall emission at 21 cm with 
differential brightness temperature of order a few mK. 
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By identifying the halos in our simulations, we were able to separate and 
compare the relative contributions of the halos and the IGM gas to the total 
signal. We find that the emission from minihalos dominates over that from the 
IGM outside minihalos, for z < 20. The minihalos contribute 70 — 80 % of the 
total emission signal at z < 16, exceeding > 85 % at z ~ 16 — 18, and balancing 
the absorption by the IGM gas at z ~ 20. In contrast, the absorption by cold 
IGM gas dominates the total signal for z > 20, with a signal close to that of 
an unperturbed IGM. 
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